Quantum probability distribution of arrival times and 
probability current density 



V. Delgado 

Departamento de Fisica Fundamental y Experimental, 
Universidad de La Laguna, 38205-La Laguna, Tenerife, Spain 
e-mail: vdelgado@ull.es 
(May 19, 1998) 



Abstract 

This paper compares the proposal made in previous papers for a quantum prob- 
abiUty distribution of the time of arrival at a certain point with the corresponding 
proposal based on the probability current density. Quantitative differences between 
the two formulations are examined analytically and numerically with the aim of 
establishing conditions under which the proposals might be tested by experiment. 
It is found that quantum regime conditions produce the biggest differences be- 
tween the formulations which are otherwise near indistinguishable. These results 
indicate that in order to discriminate conclusively among the different alternatives, 
the corresponding experimental test should be performed in the quantum regime 
and with sufficiently high resolution so as to resolve small quantum effects. 

PACS number(s): 03.65.Bz, 03.65. Ca 



I. INTRODUCTION 

The problem of incorporating the time-of-arrival concept in the theory of quantum 
measurement has remained controversial over the years, and even nowadays this ques- 
tion is open to debate [p|-[l^. In recent times this issue has acquired renewed interest 
in part due to the development of new experimental techniques for probing quantum 
systems in the time domain. For instance, by exciting an atomic system with a pulsed 
laser and measuring the subsequent flux of electrons ejected from autoionizing states, 
as a function of the time of arrival at the detector, one can gain important physical 
information which is not obtainable by probing the system in the more familiar energy 
domain jl^. On the other hand, the time domain is more related to the macroscopic 
phenomena and for this reason turns out to be particularly suitable for investigating 
quantum systems at the mesoscopic scale [pO| . 

Another related issue that has stimulated considerable theoretical effort is that 
concerning the definition and characterization of tunneling times [ ^p^ . In connection 
with this problem, Dumont and Marchioro proposed the probability current density as 



a quantum definition for the (unnormalized) probability distribution of arrival times at 
an asymptotic point behind a one-dimensional potential barrier Q. 

There exists additional motivation for trying to incorporate such a definition into 
the formalism of quantum mechanics. Indeed, the average current {J{X)) of a classical 
statistical ensemble of particles propagating in one spatial dimension along a well- 
defined direction plays the role of a probability distribution of arrival times at X. A 
simple way for translating such a result into the framework of quantum mechanics con- 
sists in invoking the Weyl-Wigner quantization rule, which provides a prescription for 
constructing a quantum operator A{X , P) corresponding to a given classical dynamical 
variable A{x,p) pS, 



A{x,p) ^ A{X,P) = j j I A{x,p)e'^^'^^-''^+<^-P^Ux dpde dr. (1) 

Furthermore, the operator so obtained has the nice property that its expectation value 
is given by the classical expression 

{A{X,P)) = llf^{x,p)A{x,p)dxdp, (2) 

with the Wigner function fy/{x,p) playing the role of a quasiprobability distribution 
function in phase space. 

The Weyl-Wigner quantization rule must be used with caution for it does not nec- 
essarily lead to the correct quantum operator. In the present context, one obtains that 
the Weyl-Wigner operator corresponding to the classical current J{X) = p/m5{x — X) 
is nothing but the usual current operator 

JiX) = ^(p\X){X\ + \X){X\p). (3) 

However, unlike the classical case, because of the fact that J{X) is not positive definite, 
its expectation value cannot be properly considered as a probability distribution of 
arrival times. It has been argued, nonetheless, that asymptotically far from a potential 
barrier the transmitted current becomes positive, and this circumstance justifies its 
interpretation as a probability distribution 0,^]. In this regard, McKinnon and Leavens 
P,p!5[ have also shown that within the framework of Bohmian mechanics it is possible 
to unambiguously define a probability distribution of the time of arrival in terms of 
the modulus of the probability current density. However, even though such a definition 
circumvents the problem mentioned above, in principle there is no justification for 
extrapolating it to the framework of standard quantum mechanics. 

A natural way for introducing time into the quantum framework as a physical 
variable consists in considering it as such already at the classical level (a fact that can 
be implemented by making a suitable canonical transformation) and then quantizing 
the corresponding formulation by using the canonical quantization method in order 
to look for the desired probability distribution in terms of the spectral decomposition of 
an appropriate self-adjoint operator. In doing so, one arrives at a time operator defined 



2 



as the operator canonically conjugate to the relevant Hamiltonian |T^,26,27|. However, 
in general, no such a self-adjoint operator exists |ll|j3|,[l2| . This is the technical reason 
that explains to a great extent the difficulty found for incorporating a time operator 
into the quantum formalism. 

A reasonable way of circumventing this problem consists in looking instead for a 
self- adjoint operator with dimensions of time not strictly conjugate to the Hamilto- 
nian. Even though there exist appreciable differences among them, the approaches of 



Kijowski H, Grot et al. [1C|, as well as the one developed in Refs. |12| and ||T3| can 



be ascribed to this category. The first two approaches are concerned with the time of 
arrival of a free particle, and its supposed range of validity includes quantum states hav- 
ing, in the momentum representation, positive- and negative- momentum components, 
while the latter is also applicable (asymptotically) in the presence of a one-dimensional 
scattering potential and its range of validity is restricted to quantum states having 
either positive- or negative-momentum contributions. In this paper we shall focus on 
this latter approach. It should be remarked, however, that within their common range 
of applicability all of them provide the same theoretical prediction for the probability 
distribution of the time of arrival at a certain point. 

Agreement with a conclusive experimental test is the ultimate requirement for estab- 
lishing the validity of any theoretical proposal. However, discriminating experimentally 
among different alternatives is not always a straightforward matter. It may happen that 
under certain experimental conditions predictions corresponding to different proposals 
become indistinguishable in practice. This is the case in the present context when con- 
sidering quantum states largely semiclassical in character. Indeed, in the semiclassical 
limit the proposal for the probability distribution of arrival times based on the 
operator approach coincides with that based on the modulus of the probability cur- 
rent density, which is the result obtained by McKinnon and Leavens within Bohmian 
mechanics [p|,[T5|. More generally, since in this limit the quantum current becomes 
necessarily positive, it follows that the predictions based on the operator approach be- 
come in fact indistinguishable from those based in general on the probability current 
density. Consequently, any experimental test performed under these particular condi- 
tions would be inconclusive. It is therefore worthwhile to investigate quantitatively to 
what extent appreciable differences among the competing proposals can be expected as 
well as to examine how such differences depend on the various controllable parameters. 
This is the main purpose of the present work . More specifically, quantitative differ- 
ences between the two formulations will be examined analytically and numerically (as 
a function of both the initial quantum state describing the particle and the parameters 
characterizing an intermediate potential barrier) with the aim of establishing conditions 
under which the proposals might be tested by experiment. To this end we shall begin 
by briefly reviewing the required formulation. 
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II. PROBABILITY DISTRIBUTION OF ARRIVAL TIMES 



Consider a quantum particle moving along the x axis toward a detector located at 
a certain asymptotic point X behind a one-dimensional scattering center V{X). In 
looking for a probability distribution of the time of arrival for such a physical system, 
we introduced in previous papers |12,13| a self-adjoint operator with dimensions of time 
T{X) defined in terms of its orthogonal spectral decomposition by 



f{X) 



\r;X) 



dTT\T;X){T;X\, 



JsgniP)£-T/h 



\P\ 



m 



\X) 



The operators sgn(P) and yl-P | are in turn given by the expressions 



sgn(P) 



dpJ \p 



dp{\p){p\- \ -p){-p\) 



(4) 
(5) 

(6) 



(7) 



where the momentum eigenstates {\p)} are assumed to be normalized as {p\p') = 5{p- 



Note that the above equations define, in fact, a one-parameter family {T(X)} of 
self-adjoint operators (labeled by the position X of the detector) which are canonically 
conjugate to the operator Ti. = sgn(P) Hq, with Hq = P"^ /2m being the energy of the 
free particle. 

Let Q{+P) [Q{—P)] represent the projector onto the subspace spanned by plane 
waves with positive [negative] momenta, 



G(±P) 



dp \ ±p){±p\ . 



(8) 



By taking advantage of the resolution of the unity 0(+P) + 0(— P) = 1, we can rewrite 
the eigenstates |t; X) (which are manifestly symmetric under time reversal) in the form 



|r; X) = e(+P) e^^«"/'* \ ^ \X) + Q{-P) g-^^o"/^ \ ^ 

\ m \ m 



(9) 



which involves the state ^ \P l/mlX) translated (freely) both forward and backward 
in time by the amount r. 

Substituting then Eq. into Eq. (Q), one obtains 



T(X) =e(+p) 



dTTjf\X, t) 



e(+p)-e(-p) 



dTTjf\X,T) 



e(-p), 

(10) 
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where the positive-definite current jj {X, r) is a straightforward quantum version (in 
the interaction picture) of the modulus of the classical current | J(X)| = \p\/m5{x — X), 

jW(X,r) = e^^"^/^ e"^^"^/'*, (11) 



jW^X)^J\£ls{X -X)J\^. (12) 
V m ^ m 

Even though T(X) is symmetric under time reversal, its restrictions to the subspaces 
spanned by either positive- or negative-momentum plane waves are not. This fact 
enables us to define a probability distribution of the time of arrival for quantum states 
belonging to either of such subspaces. To be specific, let us assume the particle under 
study to be incident from the left of the potential barrier, and let the state vector iV'in) 
[which is assumed to satisfy the identity iV'in) = 6(-P) \ ipm)] represent the incoming 
asymptote of the actual scattering state of the particle at t = 0. The mean arrival time 
at an asymptotic point X can then be defined consistently as |1^,13| 



(V'trlV'tr) {AAAr/ 

where \tptr) is the projection of the outgoing asymptote (at t = 0) onto the channel of 
transmitted particles, i.e., 

IV^tr) = e(P)IV'out) = e(P) S IVin) = / dpT{p) (plV^in) \P), (14) 

Jo 

with 5 and T(p) being, respectively, the scattering operator and the transmission co- 
efficient characterizing the potential barrier. 

It is worth noting the remarkable formal analogy between Eq. ([l^) and its corre- 
sponding classical counterpart. Indeed, the positive-definite current (V'tr I (^, ^)|V'tr) 
enters the expression for (tx) playing the role of an (unnormalized) probability distri- 
bution. We can thus define the probability distribution of the time of arrival at the 
asymptotic point X as 

Pxit) ^ j\{t;X\Ar)\'' = j{Arit)\J^+\x)\AAt)), (15) 

where T ={iptj.\il;tj.) = {iptr{t)\4^tT{t)) is the transmittance and we have written the latter 
expression in the more familiar Schrodinger picture by introducing the (Schrodinger) 
freely evolving transmitted state 

|V'tr(t)) = e-^"*/^|V^tr). (16) 

Equation (|15|) along with Eqs. ( p^ and (16) enable us to compute the desired 
probability distribution in terms of the ingoing asymptote IV'in)- It is worth, however. 
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obtaining an alternative formula in terms of the actual scattering state \ip{t = 0)). This 
can be accomplished by means of the M0ller operators, 

Q±= lim e^^*/^e-*-^o*/'', (17) 

which map the ingoing and outgoing asymptotic states onto the corresponding scatter- 
ing state 

|^(t = 0)) = = J7_|^out>. (18) 

Using these relations and introducing the projector 

r = h^e{P)nl (19) 

(which selects that part of a given state vector that will be transmitted), one finally 
obtains (Appendix A) 

Px{t) = ^mo)\'P^-ji^\x,t)Cllr\i;m = l(V'(t)|75 0_j(+)(X)AL^|V(t)), (20) 

where |V'(^)) = e~*^*'''*|^(0)) is the usual Schrodinger state vector. It is interesting 
to note that the above equation is merely the expectation value of the modulus of 
the current J^^\X) in the quantum state Cl^_V\il){t)), which, in turn, is the 

normalized outgoing asymptote corresponding to that part of \ip{t)) that is going to be 
transmitted in the future. 

In practice, whenever the actual scattering state at t = does not overlap appre- 
ciably with the potential barrier, the state vectors iV'in) and |V'(0)) become physically 
indistinguishable ||2^ and, consequently, one can legitimately use Eqs. ([T^)-(p^ with 
the substitution \tpin) — > 1^(0)). 

For our purposes it is convenient to write the expectation values of J{X) and 
J(+)(X) as 

{Ar{t)\J{X)\Ar{t)) = ^\ (r M/[l] + C.C.) , (21) 

{Ar{t)\.n^\x)\AAt)) = ^ invpinvp]) , m) 

where we have mtroduced the functional 

roo 2 

I[f]^ dpT{p) f{p) (piV'in) e-*i^*/^e^^'^/^ (23) 
Jo 

Note finally that the free case is a particular case of the above formulation with 
the M0ller operators (l± reducing to the unit operator and, consequently, P Q{P). 
Since the scattering operator can be written as S =Q,^_Cl^, it also follows that S ^i, 
and hence, by virtue of Eq. ([T^, T{p) 1 and iV'tr) iV'in) — >■ \ip{0)). With these 
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substitutions the above formulas are applicable to the study of the arrival time of a 
free particle at a point X. 



III. ANALYTICAL APPROXIMATION FOR THE EXPECTATION 

VALUE OF J(+)(X) 

In this section we are interested in obtaining analytical expressions that permit us to 
compare the proposed probability distribution of the time of arrival (^tr(i)l'^''^'*(^)lV'tr(i)) 
with the standard probability current density (V'tr(i)|'/(-'^)|V'tr(i))- To this end we shall 
restrict ourselves to a free particle characterized, at t = 0, by a minimum Gaussian wave 
packet with centroid xq, having a negligible contribution of negative-momentum com- 
ponents, and propagating with average momentum po > along the x axis toward a 
detector located at a certain position X > xq. Specifically, 



(p|V'(0))= 27r(Ap) 



-1/4 



exp 



V-Po \ _ - pxp 
2Ap J ^ h 



(24) 



where the momentum spread Ap <^ po is assumed to be sufficiently small so as to 
satisfy {p\ip{0)) ~ 0(p)(p|V'(O)) to a- good approximation. As stated above, under these 
conditions T{p) 1 and we may substitute {iptrit)) — > \ip{t)) throughout the relevant 
formulas. The integrals involved in the definition of the probability current density 
[Eqs. ( ^T|) and (p3|)] can then be easily carried out to obtain the well-known formula 



mh 



3/2 



exp 



,(Ap)2 [[X-x,)-ft\ 



1 + 4 



(25) 

As far as the probability distribution ((4){t)\J^^\X)\ijj{t)) is concerned, some additional 
simplification is needed. We shall content ourselves with an analytical approximation 



valid up to order {Ap/po) . To this end, following Grot et al. we expand the 



argument of the functional /[^y^] entering Eq. (22) as a Taylor series about the average 
momentum pQ, i.e.. 



1 + 



P-Po 



1 IP 



Po 



2po 2 V 2po 
Substitution of this expansion into leads to 



+ 



p-Po 
Po 



I[^/P\ = iVPO 



p[i]+-i[p]-^i\p']+o(^y 

2 Po 2po \Po J 



(26) 



(27) 
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The important point is that both I\p\ and I\p^] can be written in terms of /[I], yield- 
ing an expression for (V^(t)|j(+)(X)|V'(t)) that can be easily related to the probability 
current density given in Eq. (p5|). Indeed, substituting Eq. (^) in the integrand of 
/[I] and taking into account that (p|V'(0)) — for p < 0, one arrives at the Gaussian 
integral 

/ + 00 2 
^^g-5(p-A) ^ (28) 
-oo 



where 



+ (29) 



4(Ap)2 2mn 

^ Po + 2iiAp)\X-xo)/n 

l + 2i{Ap)H/mh • ^ ' 

(The factor N is not relevant to our purposes and consequently is not given here.) Note 
that Re((^) > 0, as required for the integral to converge. As is well known, the integral 
of Eq. ( |28| ) is, in fact, independent of A. Thus, by differentiating with respect to 
A one can readily show that 

i[p] = xm. (31) 

Likewise, a second differentiation with respect to A leads to 

V] = (a' + ^)/[1]- (32) 
By inserting Eq. (^) into Eq. (^I[), one obtains 

mh 

Substitution of Eqs. ( ^T| ) and (|32|) into Eq. ( |27| ) yields an expression for I[^/p\ depend- 
ing only on Inserting then the expression so obtained into Eq. (|2^ ) and using Eq. 
(p^) to eliminate in favor of the probability current (^(t)| J(X)|^(t)), we arrive 

at 



(^(t)|j(+)(X)|^(t)) 



+ n'^' + + ° (^)l W'(«)|i(X)l*(*)>^(34) 



Of course, for this cumbersome expression to be useful some simplification is still re- 
quired. Let to = (^ ~ XQ)m/pQ be the classical time of arrival at the detector located 
at X. Restricting ourselves to particles arriving in the time interval [0,2to], we have 
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t < 2to 



{Apft 
mh 



<P 



(35) 



Po J 

where the parameter p denotes the distance between the centroid of the wave packet 
at t = and the detector's position, in units of the spatial spread Ax = h/2^p, i.e., 
{X — xq) = yoAx. The important point is that when p ~ 0(1), the term 2{Ap)'^t / mh 
involved in the expressions of both 1/6 and A [Eqs. (|29|) and (^0|)] becomes of order 
(Ap/po). This fact enables us to approximate the various terms contributing to Eq. 
(p^) as Taylor's expansions up to terms of order (Ap/po)^- For instance, under the 
assumptions just stated we would have 



Re(A) = Po 



1 + 4 



{Ap)^iX-xo)t 
mpoti'^ 



\Po J 



(36) 



with similar expansions for the rest of the terms entering the above expression for 
(V'(t)|J^+H-^)|V'(i)>- The substitution of all of these expansions into Eq. (|3^ ) leads, 
after a rather lengthy calculation, to the final result 



{m\j^^\x)\m) 



l + Ao- Ait + A2t" + ( 



m)\j{x)\m), (37) 



where the coefficients Aq, Ai, and A2 are defined as 



(38) 

2 \Po J fi Pq 

mpoh 

M . 2<^. (40) 

Equation ( |37|) constitutes the main result of this section. When used along with Eq. 
(p^), it enables us to obtain analytically the probability distribution (i/^(t)| {X)\^l^(t)) 
up to order {Ap/po)^. When the detector is located at a position X initially separated 
from the centroid xq a distance of order Ax, its range of applicability extends over 
the time interval [0, 2to]. However, the validity of Eq. (^) is not restricted to this 
particular configuration. Indeed, it is not difficult to see that for any detector location 
X > Xq (i.e., any p > 0), Eq. (^7|) remains true within the interval [0,crto], with 
a ~ min {(2//)), (2/p)^}. For greater times the contribution from those terms neglected 
in the various expansions increases, so that they can no longer be considered to be of 
order (Ap/po)^ and Eq. (^) might fail as an approximation valid up to order (Ap/po)^- 

As is apparent from the above formulas, the difference between the expectation 
values of J{X) and J^'^\X) turns out to be of order {Ap/pof', so that, in the limit 
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(Ap/po) 0, both of them coincide, a result that is in agreement with Ref. |1C] 
and confirms the asymptotic analysis of Ref. [^]. A comparison between these two 
quantities is, however, most conveniently done in terms of their relative difference. 



IV. RELATIVE DIFFERENCE 

Given a Hamiltonian H = Hq + V{X) and a state vector having no contri- 

bution of negative- momentum components, and not overlapping appreciably with the 
potential barrier [so that one may legitimately substitute IV'in) — * lV'(0))]i the relative 
difference A between the probability distribution {il;tr{t)\J^^\X)\'ilJtr{t)) and the cor- 
responding probability current density (V'tr(OI'^(^)lV'tr(i)) can be written as [see Eq. 
(0)] 



^ ^ ^ _ {m\s^QiP)Mx,t)eiP)s\m) (4^) 
(V^(o)|5te(p)j;+)(x,t)G(p)5|^(o))' 

This quantity, which constitutes the basis for our subsequent analysis, is (for any t) a 
state functional that enables us to quantify the differences we are interested in, as a 
function of the initial state \ip{0)). To this end we shall consider an initial state vector 
of the form 

1^(0)) = a(/3|Vi) + |V2», (42) 
where a is the normalization constant, 

a = (/32 + 2^Re(V'i|V'2) + (43) 

and P is an arbitrary real coefficient. On the other hand, for computational simplicity 
we shall choose {p\^j) (j = 1, 2) to be minimum Gaussian wave packets centered at xq, 
with momentum spread Ap and average momentum pj, respectively. 



{pm= 27r(Ap) 



-1/4 

exp 



2Ap J * n 



(44) 



Furthermore, we take P2 ^ Pi > and Ap <^ pi in order to guarantee that \ip{^)) has 
no appreciable contribution of negative-momentum components. 

The interest in choosing \ip{Q)) this way comes from the fact that by varying the 
continuous parameters (i, pi, p2, Ap, and xq we can easily explore different regions 
of the Hilbert space. Such an analysis will be the aim of the next section. For the 
time being, consider a free particle described at t = by the state vector defined by 
Eqs. (^2|)~(|4^) withp2 Pi = Po- In this particular case, |V'(0)) — > IV'i) and the initial 
state reduces to the simple minimum Gaussian wave packet considered in the preceding 
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section [Eq. (24)]. Under these circumstances, an analytical expression can be derived 



for A. By substituting Eq. (§3) into Eq. @, one finds (for < t < ato) 

A = A.(t-*o)^-im%of^)'. (45) 



2 V Po / \Po J 

with A2 given in Eq. (|40|). Therefore, as a function of t, the relative difference A is given 
by a parabola which reaches its minimum at the classical arrival time to = {X—xo)m/pQ. 
From Eq. (]45| ) it follows that A (to) < 0, and consequently 

(V^(to)|jW(X)|^(to)) < mo)\JiX)mto)), (46) 

so that the probability of arriving at to as predicted by the current J{X) is always 
greater than that predicted by the modulus of the current J^~^\X). More generally, it 
can be readily seen that for any t > within the symmetric interval [to — mAx/po, to + 
mAx/po] about to, it holds that 

(v(t)|j(+)(x)i^(t)) < {m\j{x)\m), (47) 

where the equality is satisfied at the boundaries of the interval [up to terms of order 



V. QUANTITATIVE ANALYSIS 
A. Free particle 

We begin the analysis of the relative difference A between the probability distribu- 
tion {'il^t,-{t)\J^~^\X)\\l)t,-{t)) and the probability current density (V'tr(i)|</(-'^)|V'tr(0) by 
considering a free particle described at t = by the minimum Gaussian wave packet 



defined in Eq. (24) [which, as already said, is nothing but a particular case of the 
state vector previously introduced in Eqs. (^2[)-(^)]. To be specific, we shall restrict 
our investigation to an electron with average momentum pQ = 0.5 a.u., and place the 
detector at X = xq + 3Ax, where Ax = h/2Ap is the spatial spread of IV'(O)). 

Figures 1(a) and 1(b) show the relative difference A, as a function of the time of 
arrival, for several values of Ap ranging over two orders of magnitude. The detection 
interval has been chosen to be symmetric about the classical arrival time to, i.e., t G 
[0, 2to]. Moreover, for comparison purposes the final instant of time has been normalized 
to unity in all cases (i.e., t ^ tn = t/2to), so that to lies always in the middle of the 
detection interval considered. It is worth remarking that the probability current density 
turns out to be positive in all of the cases studied. 

Besides the results obtained from a direct numerical integration of the corresponding 
formulas, Fig. 1(b) also shows the behavior of A as predicted by Eq. (p5|). (Only the 
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theoretical values corresponding to Ap = 0.01 a.u. have been explicitly plotted [circles 
in Fig. 1(b)] since for smaller Ap the agreement is even better.) As is apparent 
from this figure, within the range of validity of the theoretical prediction (that is, for 
(Ap/po) <C 1 and t £ [0, ato]) the agreement turns out to be excellent (as it should be). 
In fact, for any momentum uncertainty Ap < 0.01 a.u., the relative difference between 
(j(+)(X))(t) and {J{X)){t), considered as a function of t, behaves clS Si parabola which 
cuts the horizontal axis at = 1/3 and = 2/3, reaching its minimum value at 
the classical time = 1/2. Furthermore, while within the interval [1/3,2/3] the 
expectation value of J{X) always dominates over the expectation value of J^^\X), 
just the opposite occurs outside such an interval. 

From Figs. 1(a) and 1(b) the rapid decrease of A with the momentum spread is 
also apparent, which is a direct consequence of the fact that for Ap/pQ sufficiently 
small, A is of order {Ap/pof' . In this regard, note that for Ap < 0.01 a.u. the relative 
difference is already less than 0.2% over all of the detection interval considered, and that 
figure would be even smaller if one focused attention on a time interval more localized 
about the most probable time of arrival to- This fact would render any attempt to 
discriminate between {J^^\X)){t) and {J{X)){t) almost impossible in practice. Such 
a negative conclusion should not be extrapolated, however. A minimum Gaussian wave 
packet represents a very special type of quantum state, for it exhibits the lowest possible 
uncertainty product AxAp and consequently is expected to be largely semiclassical in 
character. On the other hand, as can be inferred from an asymptotic analysis [13], in 
the semiclassical limit h ^ it holds that 



so that the small value found for A in the cases considered above is not surprising. 

A simple way for generating initial states having a more genuine quantum charac- 
ter consists in allowing pi to be different from p2 in Eqs. (^)-(p^). In fact, a marked 
interference pattern can be induced in the probability current density (as well as in the 
corresponding position probability density) by simply increasing the distance (in mo- 
mentum space) between pi and p2 while keeping unchanged the remaining parameters. 

Figures 2 and 3 show the results obtained by taking (3 = 2, p2 = 0.5 a.u., Ap = 0.01 
a.u., and allowing pi to vary between 0.4 and 0.2 a.u., respectively. The detector's 
position has been chosen as before (more precisely, we have taken X — xq = 2>h/2Ap) 
and the detection interval (now expressed in atomic units) has been chosen in such a 
way that its final time tf satisfies 

m,)\j^+\xmt,)) « (V'(o)ij(+)(x)iV'(o)). (49) 

In Figs. 2(b) and 3(b) we have plotted the expectation values of both j'^'^\X) and 
J{X) corresponding, respectively, to pi = 0.4 a.u. and pi = 0.2 a.u. These curves 
can be considered as being obtained by means of a continuous deformation (induced by 
varying pi) starting from the initial Gaussian profile corresponding to pi = P2 = 0.5 a.u. 
Since a marked interference pattern is a hallmark of quantum behavior, it is evident 
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that by decreasing pi we are probing domains of the Hilbert space with increasing 
quantum character. As is apparent from these figures, while for pi = 0.4 a.u. the 
probabihty current density remains positive over ah of the time interval considered, 
the same does not occur for pi = 0.2 a.u., and in this case {ip{t)\J {X)\ip(t)) attains 
negative values in the neighborhood of t = 40 a.u. and t = 400 a.u. 

The analysis of the corresponding relative differences, plotted in Figs. 2(a) and 
3(a), reveals that now A exhibits a series of maxima (clearly related to the interference 
pattern of the probability current density) where it can take values of order 10% (in the 
first case) or 100% (in the second one). [In fact, in this latter case the maxima which 
have been truncated in the figure correspond, respectively, to A(39.1 a.u.) = 31.96 
and A(398.05 a.u.) = 1011.58.] Therefore, the probability distributions of the time of 
arrival as predicted by {J^+\X)){t) or {J{X)){t) can be quite different in these cases. 
A comparison with the bound 0.2% obtained previously for pi = P2 = 0.5 a.u. [Fig. 
1(b)] reflects the fact that the choice of the initial state plays an essential role in both 
the magnitude and the behavior of the relative discrepancy between (j(+)(X))(t) and 
{JiX)){t). 

An interesting limiting situation, where quantum effects dominate the behavior 
of the probability current density, can be achieved by taking Ap — > and P2/P1 3> 



/3 ^ 1 in Eqs. (|4^)-(44). In this particular case, the initial state becomes a purely 
quantum state with no classical analog, consisting of the coherent superposition of two 
macroscopically distinguishable states in momentum space. Under these circumstances, 
the main contribution to {•ip{t)\J {X)\'ip{t)) comes from the interference terms and one 
obtains, to leading order as Ap ^ flc 



2 /27r 

{iP{t)\J{X)\7p{t)) ^ ^^(3a^Ap P2COS plt/2hm-p2{X -xo)/n +0 {Apf . (50) 
mh L J L J 

This case is illustrated in Figs. 4(a) and 4(b), where we have specifically taken (3 = 100, 
P2 = 1 a.u., pi = 4x 10~^ a.u., Ap = 5 x 10~^ a.u., and X — xq = 2>h/2Ap. Also plotted 
in Fig. 4(b) is the modulus of the probability current density, \ {ilj{t)\J {X)\'ilj{t))\, which 
is the proposal by McKinnon and Leavens ||8| for the probability distribution of the time 
of arrival. 

Prom a comparison between Figs. 4(a) and 4(b) it can be seen that the relative 
difference A between {J{X)){t) and {J^'^\X)){t) reaches a local maximum exactly 
when {J{X)){t) reaches a local minimum, and in this case A S> 1, so that 

{j{x)m \{j{x)m\ 

(j(+)(x))(t) (J(+)(x))(t)- 

Consequently, the main effect produced by the replacement {J{X)){t) — > \{J{X)){t)\ 
upon the corresponding relative difference would consist in a change of sign. A similar 
conclusion could be inferred from Fig. 3(a) for those instants of time where the prob- 
ability current density becomes negative. The important consequence is that, from 
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a quantitative setting, by replacing the probability current density by its modulus 
one in general does not achieve a better agreement with the probability distribution 

(J(+)(X))(0. 



B. Potential barrier 

The presence of a potential barrier might, in principle, induce some additional 
discrepancy between {J{X))(t) and {J^~^\X))(t). In order to examine whether this 
is the case, we shall consider next a quantum particle propagating toward a detector 
located at a certain asymptotic point X behind a one-dimensional potential barrier. 
It turns out to be most convenient restricting to initial states as simple as possible 
since in these cases the effect produced specifically by the barrier can be more easily 
identified. Accordingly, we shall consider an electron characterized at t = by the 
Gaussian state given in Eq. (p4[), with pQ = 0.5 a.u. The potential barrier is assumed 
to occupy the segment [0, d] of the x axis, and its height Vq has been chosen in all cases 
so that pb = \/2mVo = 0.8 a.u. 

Our primary interest consists in investigating the behavior of the relative difference 
A as a function of both the barrier's width d and the momentum spread Ap. To this 
end we have selected the centroid xq of the initial wave packets in such a way that 
{x\ip{0)) does not overlap appreciably (in comparison with the transmittance T) with 
the interaction center. Specifically, xq has been implicitly defined by 

roc 

/ dx\{x\ip{0))\'^^lO-^J (52) 
Jo 

and, consequently, is a function of both Ap and d. On the other hand, in order to 
guarantee the applicability of the formalism (which requires the particle to be asymp- 
totically free), the detector has been assumed to be switched on at a certain instant ti 
satisfying the condition that the probability of finding the particle within the interac- 
tion region is already negligible. More precisely, we have defined implicitly ti by the 
condition 

r dx\{x\Arim^ ^W^^T. (53) 

The detector has also been assumed to be located at a certain position X behind the 
interaction center sufficiently far from the barrier's edge so as to satisfy the condition 
that when it is switched on (at t = ti) the probability of finding (in the absence of 
detector) the transmitted particle within the region x > X is still negligible. That is, 
X has been obtained from the condition 

POO 

/ dx\{x\Ar{ti))\^ (54) 
Jx 
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Table 1: Parameters corresponding to Ap = 0.01 a.u. All values in atomic units. 



d 


Xq 


ti 


if 


X 


2 


-201.8 


785 


1550 


379.0 


4 


-228.0 


839 


1600 


382.0 


8 


-275.4 


933 


1700 


386.5 


12 


-316.6 


1014 


1800 


391.0 



Finally, it has been assumed that the detector is switched off at a certain instant t{ 
satisfying 

{AAti)\j'^^HX)\AAti)) ^ {Ar{ti)\J^^HX)\Ar{ti)). (55) 



Combining Eqs. (|5^) and ( |54| ) we see that at t = t^, i.e., when the detector is 
switched on, the transmitted particle can be found in the region between the barrier's 
edge and the detector's position with a probability of 99.8%, a fact that guarantees the 
applicability of the formulation developed in the preceding sections. 

Figure 5 shows the relative difference A obtained by taking Ap = 0.01 a.u. and al- 
lowing the barrier width to vary between d = 2 a.u. and d = 12 a.u. The corresponding 



values for the parameters xq, t[, tf, and X [obtained numerically from Eqs. (52)-(55)] 
are given in Table ||. For comparison purposes, the time interval [ti, if] has been renor- 
malized to the interval [0, 1] by means of the mapping t ^ = {t — ti)/(tf — t]). As is 
apparent from Fig. 5, the potential barrier seems to produce no special effect on A. In 
fact, even though the transmitted state |V'tr(i)) depends to a great extent on the barrier 
width d, the relative difference between the corresponding expectation values of J{X) 
and J^^\X) exhibits no appreciable dependence on this parameter. Furthermore, A 
remains very small (of order 0.25%) in all cases. 

Things are quite different, however, when the momentum spread increases. For 
Ap = 0.1 a.u. (Figs. 6 and 7 and Table §), the relative difference not only takes greater 
values, but also exhibits a clear dependence on the barrier width. From Fig. 6(a) we see 
that for d = 8 a.u., A manifests an oscillatory behavior which [from a comparison with 
Fig. 6(b)] can be straightforwardly related to the emergence of an incipient interference 
pattern in the corresponding probability current density. The appearance of such an 
interference pattern can in turn be traced back to the interference dynamically induced 
between tunneling and over-the-barrier contributions in the corresponding transmitted 



wave packet |3C]. 



A similar behavior associated, however, with greater relative differences can be ap- 
preciated in Figs. 7(a) and 7(b), which show the results corresponding to a barrier 
width d = 10 a.u. This particular case is interesting for still another reason. As 
is shown in the inset of Fig. 7(b), the probability current density takes in this case 
negative values in the neighborhood of the arrival time t = 373 a.u. (though, admit- 
tedly, very small ones). This fact demonstrates that even for initial Gaussian wave 
packets having no appreciable contribution of negative-momentum components, the 
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Table 2: Parameters corresponding to Ap = 0.1 a.u. All values in atomic units. 



d 




ti 


ti 


X 


2 


-20.15 


145.7 


530.0 


112.9 


4 


-22.48 


137.0 


470.0 


108.5 


8 


-25.84 


141.1 


390.0 


117.3 


10 


-28.30 


254.0 


573.0 


229.4 



probability current density can take negative values at an asymptotic point X behind 
a one-dimensional potential barrier. This sole reason is sufficient to invalidate (even 
under such special circumstances) its interpretation as a probability distribution of 
arrival times. Of course, there still exists the possibility of considering instead the 
modulus of the probability current density. Were we to plot the relative difference 
between {J^^\X))(t) and \ {J{X)){t)\, the only relevant change in Fig. 7(a) would be 
the transformation of the first maximum A(t = 373.3 a.u.) = 3.92 into a minimum 
(3.92 —1.92). Consequently, a considerable discrepancy would still survive. 



VI. CONCLUSION 

In the present work we have compared the proposal made in previous papers |12,13| 



for a quantum probability distribution of the time of arrival at a certain point with that 
based on the probability current density (or alternatively on its modulus), with the aim 
of establishing conditions under which the proposals might be tested by experiment. 
To this end, we began by obtaining (under certain particular conditions) an analytical 
approximation for the expectation value of J^~^\X) valid up to order (Ap/po)^, and 
we have performed a quantitative analysis of the corresponding relative differences as 
a function of the initial state of the particle (both in the case of free evolution and in 
the presence of an intermediate potential barrier). 

We have found that quantum regime conditions produce the biggest differences 
between the formulations which are otherwise near indistinguishable. In fact, in the 
semiclassical regime, for electrons in quantum states with a well-defined momentum 
Po ^ Ap (thus having a clear classical analog), the relative discrepancy can be typi- 
cally of order 0.2%, and this figure would be even smaller if one restricted attention to 
arrival times having an appreciable probability. Therefore, in this regime the probabil- 



ity distribution proposed in Refs. |12,13| becomes indistinguishable in practice from the 
corresponding probability current density. Important discrepancies only occur in the 
purely quantum regime, for states of genuinely quantum character, having no classical 
analog. Indeed, the appearance of important relative differences can be straightfor- 
wardly related to the existence or emergence of a marked interference pattern in the 
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probability current density (or alternatively in the position probability density), irre- 
spective of the fact that such quantum interference was already present in the initial 
state or was generated dynamically through the temporal evolution of the system (as 
can be the case in the presence of a scattering potential). Furthermore, a closer anal- 
ysis reveals that such quantum effects are mainly localized about arrival times having 
a negligible probability and/or occur over short time scales in comparison with the 
relevant time interval. These results indicate that in order to discriminate conclusively 
among the different alternatives, the corresponding experimental test should be per- 
formed in the quantum regime and with sufficiently high resolution as to resolve small 
quantum effects. Hopefully, the recent advances in the development of experimental 
techniques as well as in the preparation and manipulation of atomic systems will make 
such experiments feasible. 
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APPENDIX A 

Let us introduce the state \p—) = r2-|p), which is the solution of the Lippmann- 
Schwinger equation corresponding to an outgoing plane wave i.e., 

\p-) = \p) + (y/2m - iO - V{X)\p). (56) 

In terms of these states one can define the projector 



P = n^Q{P)nl = l dp\p-){p-\, (57) 



oo 







which selects that part of a given state vector that will be finally transmitted. Taking 
into account that = 1 and S =Cl^_Cl^, it follows from Eq. (|57| ) that 

nlrh+ = e{p)s. (58) 



By using the intertwining relations for the M0ller operators |2£], 

n}^Hn± = Ho, (59) 
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as well as Eq. ([5^), it can be readily shown that 

Taking advantage of this relationship and recalling that IV'tr) = &{P)S {tpm) , we can 
write 

\Ar{t)) = e-^^«*/^ e{P)s iV'in) = nlv \m), (61) 

where we have used that e~*^*/^ A+l'i/'m) = e~*^*/^|'i/'(0)) = \i^{t)). In particular, for 
t = 0, Eq. ( |6l| ) reduces to 

\Ar) = ^lrmO)). (62) 

Substitution of Eq. (|5) [or alternatively Eq. (|6|)] into Eq. ([|) leads to Eq. @. 
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Figure captions 



Fig. 1. Relative difference A as a function of the (normalized) time of arrival = 
t/2to [where to, which depends on Ap, is the classical arrival time: to = {X — xo)m/pQ = 
3mfi/2pQAp]. The circles superimposed in the case Ap = 0.01 a.u. correspond to the 
theoretical prediction of Eq. (^). 

Fig. 2. (a) Relative difference A as a function of the time of arrival t. (b) Expec- 
tation values (indistinguishable on the present scale) of J^^\X) and J{X) (in a.u.) as 
a function of t. 

Fig. 3. (a) Relative difference A as a function of the time of arrival t. (b) Expec- 
tation values of J^^\X) and J{X) (in a.u.) as a function of t. 

Fig. 4. (a) Relative difference A between (j(+)(X))(t) and {J{X)){t) function 
of the time of arrival t. (b) {J^+\X)){t), {J{X)){t), and \{J{X)){t)\ (in lO'^ a.u.) as 
a function of t. 

Fig. 5. Relative difference A as a function of the (normalized) time of arrival 
tn = (t — ti)/{t{ — ti), for Ap = 0.01 a.u. and barrier widths d = 2, 4, 8, and 12 a.u. 

Fig. 6. (a) Relative difference A as a function of the (normalized) time of arrival 
in = (t — ti)/iti — ti), for Ap = 0.1 a.u. and d = 2, A, and 8 a.u. (b) Expectation 
values (indistinguishable on the present scale) of J^^\X) and J{X) (in 10~^ a.u.) as 
a function of t^, for Ap = 0.1 a.u. and d = 8 a.u. 

Fig. 7. (a) Relative difference A as a function of the time of arrival t, for Ap = 0.1 
a.u. and d = W a.u. (b) Expectation values (indistinguishable on the present scale) of 
J(+)(X) and J{X) (in 10"^ a.u.) as a function of t, for Ap = 0.1 a.u. and d = W a.u. 
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